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Abstract 

We develop a set of heuristic arguments to explain several results on 
planar Poisson-Voronoi tessellations that were derived earlier at the 
cost of considerable mathematical effort. The results concern Voronoi 
cells having a large number n of sides. The arguments start from an 
entropy balance applied to the arrangement of n neighbors around 
a central cell. It is followed by a simplified evaluation of the phase 
space integral for the probability p„ that an arbitrary cell be n-sided. 
The limitations of the arguments are indicated. As a new application 
we calculate the expected number of Gabriel (or full) neighbors of an 
n-sided cell in the large-n limit. 
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Figure 1: Voronoi tessellation of a set of point particles. Voronoi cells are convex 
and their edges join at trivalent vertices. 

1 Introduction 

Planar cellular structures occur in a wide variety of natural systems. The 
examples most quoted are certain biological tissues and soap froths confined 
between parallel plates. In addition, cellular structures are employed as a 
tool of analysis in a great diversity of problems throughout the sciences and 
beyond. Many references may be found e.g. in Okabe et al. [1] and in Rivier 

The Voronoi tessellation is one of the simplest mathematical models of 
a cellular structure. In two dimensions it is obtained by distributing a set of 
point-like "seeds" in the plane and then performing the Voronoi construction 
as in the example of Fig.[TJ the plane is partitioned into cells such that each 
point of the plane is in the cell of the seed to which it is closest. Voronoi 
cells are convex and their edges join at trivalent vertices. In the special case 
that the seed positions are drawn randomly from a uniform distribution, one 
speaks of a Pozsson-Voronoi tessellation. 

The statistical properties of the Poisson- Voronoi tessellation were studied 
by Meijering [3] as early as 1953. Among the quantities of greatest interest is 
the sidedness probability pn, defined as the fraction of cells that are ?7,-sided 
[n = 3,4,5, . . .). As n increases, pn passes through a maximum at n = 6 
and falls off to zero very rapidly for n ^ 10. The statistical properties of 
Voronoi cells range from fairly easy to very hard to determine, depending on 
the quantity of interest. Many known properties are discussed and/or listed 
in Ref. P (see, e.g., Ch. 5, Table 5.5.1). 
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It so happens that the sidedness probabihty pn is very hard to calculate. 
Although it may readily be represented as a 2n- dimensional integral, the 
variables of integration are coupled in such a way that this is a true many- 
particle problem. As a result, no simple analytic result for p„ is known for any 
n. We showed earlier, however, that progress can be made if one considers 
the limit of large n. In that limit Pn is given by [H U] 



where C = 0.344 347... Alternatively Eq. (11. ip may be written 

logp„ = -2nlogn + nlog(27rV) - ^ logn - i log (2VC~2) +o(l). (1.2) 

The interest of this asymptotic result, obtained by lengthy and technical 
mathematics, goes much beyond the formula itself. Its derivation has brought 
along a full analysis of the configurational statistics of the n-sided cell, which 
in turn has opened the door to further developments in various directions 
that we will briefly mention now. 

(i) First, let m„ be the average sidedness of the neighbor of an n-sided 
cell. In the wake of the work of Ref. [5] it was shown [7] that this sidedness 
pair correlation has the large-n expansion m„ = 4 + 37r2n~2 + . . .. This 
demonstrated that Aboav's celebrated law [8], usually written as nm„ = 
an + b, is in fact a linear approximation useful for, but limited to, small n 
values. 

(ii) Secondly, the development of a new Monte Carlo method [9] made it 
possible to simulate n-sided cells and to numerically determine pn with four- 
digit precision for all finite values of n. This method was used to generate 
snapshots of extremely large cells with sidednesses as high as = 1600 [TO]. 

(iii) Thirdly, the analytic methods developed for Voronoi cells in Refs. [H 
[5] proved to be of wider use. In Ref. [Hj they were applied to line tessellations 
of the plane and linked up with work in mathematics by Hug and Schneider 
[T2] . In Ref. [13] they were brought to bear on Sylvester's question [T3]: 
what is the probability p* that n points chosen randomly from a uniform 
distribution in a disk, are the vertices of an n-sided convex polygon? This 
question leads to a connection between on the one hand our studies of many- 
sided Voronoi cells, and on the other hand work on random points in convex 
position by Barany [15j, Barany et al. [TE], and Calka and Schreiber [17], as 
well as work in statistical physics on extremal statistics of random walks by 
Gyorgyi et al. [18j and Majumdar and Comtet [T9| [20]. 

Because of all these ramifications we are justified to ask which of the re- 
sults on many-sided Voronoi cells can be understood without the full mathe- 
matical apparatus that was necessary in Ref. [5] and in most of the following 




n 



oo. 



(1.1) 
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articles. This paper answers that question. It is a successor to Ref. [5], of 
which we will reproduce part of the results by new, elementary but heuristic 
arguments. We also indicate the limitations of these arguments. Finally, the 
present paper will lay the basis for a study of Voronoi cells in higher spatial 
dimensions |21j . 

Our presentation is as follows. A preliminary observation in section [2] has 
to do with the difference between the statistics of large and of many-sided 
cells. We then estimate Pn in two successive stages. First, in section [3l and 
considering the many-sided cell as circular, we present heuristic arguments 
leading to the simplest of all possible estimates. Secondly, in section HI and 
building on what we have learned, we perform part of the mathematics of 
Ref. [5J in a strongly simplified way. This leads to a much improved estimate 
of pn- 

In section [5] we consider an effect not captured by the preceding argu- 
ments, namely the "elastic" deformations of the cell perimeter from circular- 
ity. The key to their analysis is the so-called "random acceleration process," 
also known from the theory of Brownian motion. We derive the scaling with n 
of the elastic deformations; it is not possible, however, to obtain heuristically 
their precise contribution to p„. 

In section 15.31 finally, we consider the average number of Gabriel neigh- 
bors of an n-sided cell and show that for large n it scales as ~ raa. Adjacent 
Voronoi cells are called Gabriel neighbors when the line segment that con- 
nects their seeds does not intersect a third cell. Gabriel neighbors have been 
studied in mathematics [22l [23] and also play an important role in certain 
pattern recognition algorithms where "decision boundaries" must be con- 
structed [211 [25] . Section [6] is our conclusion. 

2 A preliminary observation 

We consider a Poisson- Voronoi tessellation with seed density p; a measure 
for the typical interparticle distance is therefore £ip = p~2. For all questions 
of interest this density may be scaled to unity, but we will keep it as a check 
on dimensionality. We select an arbitrary seed, let its position be the origin, 
and are interested in its Voronoi cell, which we will call the "central cell." 
This cell has the same statistical properties as all others. Our focus is, first 
of all, on its sidedness probability 

2.1 Large cells 

It is useful to begin by pointing out a distinction. When the central cell is 
large, in the sense of having a large area, it will also typically have many 
sides, and inversely. Nevertheless, the statistical subensemble of central cells 
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Figure 2: Black dots: seeds randomly distributed in the plane with an a priori 
uniform density p = 1. The picture shows a very low-probability configuration 
consisting of a seed in the origin surrounded by a disk of radius 2r void of other 
seeds. The Voronoi cell of the central seed will be a convex n^-sided polygon, 
where rir is a random number of order ra . The cell perimeter will roughly coincide 
with the circle of radius r shown in the figure. 

with a prescribed area A is inequivalent to the one where it has a prescribed 
number n of sides, even in the limit of large A and n. We will now discuss 
why. 

Let us first ask how it can happen that the central cell (say an approxi- 
mately circular one |26]) has a large area A = vrr^. As shown in Fig. [21 this 
occurs when a disk of radius 2r around the central seed is free of other seeds. 
Now let us inquire about the sidedness of this large cell, or, equivalently, 
its number of first-neighbor cells. One might think naively that in the limit 
of large r this number can be estimated as the product of the seed density 
p and the area of an annulus of width i^p along the perimeter of the empty 
disk. That would give n,. ~ p x £ip x Airr ~ p^r, where the sign ~ denotes 
asymptotic proportionality as r gets large. The Voronoi construction, how- 
ever, does not confirm this linearity between Ur and r. Thanks to work by 
Calka [2H] and Calka and Schreiber [T7| we know that in fact the number of 

1 2 

first neighbors increases only as rir ~ (par) a and that they are essentially 

2 1 

located in a circular annulus of width ~ p~3r~3 along the perimeter of the 
empty disk. The cause is a screening effect that we will study in detail in 
section 13. 1[ 
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2.2 Constructing a many-sided cell 

Obtaining an n-sided cell by the construction of Fig. [2] would require having 
an empty disk of radius 2r such that = n, and hence would require r ~ 
Tn = p~2n2 . The probability for this to happen is exp(— vrpr^) = exp(— est x 
n^), which with respect to an arbitrary seed distribution corresponds to an 
entropy loss AS equal to [29] 

AS r^^ n-^oo. (2.1) 

It then becomes clear that perhaps there is a way of creating an n-sided 
central cell at lower entropy cost. If the annulus with the n first neighbors 
is contracted homothetically, the central cell remains n-sided. The reduction 
of the annular surface costs entropy but there is a compensation due to the 
extra space that becomes available to the seeds outside the annulus. In the 
next section we will obtain the maximum-entropy arrangement by balancing 
these two effects against each other in what may be called an entropy-versus- 
entropy argument. 

3 Heuristic estimate for Pn 

In order to heuristically estimate the sidedness probability p„ we proceed in 
two stages: a description of the screening effect in section [STT] and an entropy 
balance in section 13. 2[ 

3.1 Screening 

Fig. [3] shows the central seed in O together with five of its first-neighbor 
seeds located at Pi, . . . , P5. By constructing in the "mid-points" Ri, . . . , R5 
the perpendicular bisectors of OPi, . . . , OP5 one obtains five segments of the 
perimeter of the central cell. The positions of the vertices S2, ■ ■ ■ , S5 then 
follow as indicated in Fig. HI We will refer to the line interval SmSm+i as 
the mth perimeter segment (with cyclic boundary conditions, Sn+i = Si). 
Obviously, integrating over the n nearest-neighbor positions Pi, . . . , P„ is the 
same thing, apart from trivial factors 2, as integrating over the n mid-point 
positions Pi, ... , P„. 

We ask now the following question: what is the area available to an 
arbitrary seed if all other seed positions as well as the number n of perimeter 
segments are kept fixed? Specifically, let us study in Fig.|l]how P3 can move 
around when Pi, P2, P4, P5 are fixed. Both the polar angle and the length of 
the vector OP3 may vary, but S'3S'4 should stay perpendicular to this vector. 
The domain of variation of P3 is, clearly, determined by the condition that 
the intersection points 6*3 and 6*4 remain confined to the immobile intervals 
S2T and S^T, respectively. Some reflection then shows that P3 is restricted 
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Figure 3: A central seed located in the origin O and five of its neighboring seeds at 
positions Pi, . . . , P5. The "mid-points" . . . , -R5 are the centers of the intervals 
OPi, . . . , OP5. Heavy solid line: part of the perimeter of the central Voronoi cell. 
Heavy dashed lines: boundaries between the first-neighbor Voronoi cells. 




Figure 4: Heavy solid line: same part of the perimeter of the central Voronoi cell 
as shown in Fig.[3j The vertices are denoted S2, ■ ■ ■ , S^; point T is a virtual vertex: 
it would be there if the perimeter segment S^S/^ were absent. The shaded region 
is the area available to mid-point R3 when all other mid-points remain fixed and 
the number of cell sides is not allowed to change. This region is bounded by three 
circular arcs; the circles pass through the origin and have their centers in S'2, T' , 
and halfway the intervals 05*2, OT, and OS^, respectively. 
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Figure 5: The three circles whose arcs were shown in Fig.UJ are represented here 
in full and for much larger n. Their diameters are approximately equal to R and 
their centers form angles or order with the origin. In the limit of large n the 
shaded area becomes triangular with two acute angles of order n~^. 

to the shaded region in Fig.|H This region is bounded by three arcs of circles 
centered at points 5*2, T', and S'^ which are located halfway between the 
origin and the corresponding vertices. The exact surface area of the shaded 
region may be expressed in terms of the mid-point coordinates, but we will 
not need that much detail. 

In the limit of large n, shown in Fig.[5l the angles S2OT' and S'^^OT' will 
typically be of order n~^. The three circles have been represented in full in 
this figure. Let R be their typical diameter (which is also the radius of the 
central cell). One easily shows that for large n the shaded area tends towards 
a triangle with a long side of order Rn~^ and two acute angles of order n~^. 
It follows that the triangle height w scales as 

w ~ Rn^^, n —y 00. (3.1) 

This equation expresses the screening effect: the larger n, the narrower the 
annulus containing the n points. We will make essential use of relation (13. ip 
in the subsection 13. 2[ It is worth noting the following consequence. The 
shaded area in Fig.H] shows that when the polar angle of R^ approaches the 
one of i?2 or R4, the available phase space goes to zero linearly in the angle 
difference. This means that the screening leads to an effective repulsion 
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Figure 6: Given an annulus of outer radius 2R and width w surrounding a central 
seed, we ask for the probabihty that out of N other seeds in an area L^, there 
be none in the inner disk of radius 2R — w and exactly n inside the annulus, the 
remaining ones occupying the outer region. 

between the mid-points, an effect that will be made more quantitative in 
section HI 



3.2 Entropy balance 

The considerations above suggest the next question. We consider a central 
seed in the origin and others in an area L x L. At a suitable point we 
will let N,L ^ oo with p = N/L^ fixed. We now ask what the probability 
is for the seed configuration to be constrained as illustrated in Fig.O of all 
seeds, n are located in an annulus of outer radius 2R and width w, and the 
remaining ones are all outside of this annulus. This probability will be our 
approximation for the sidedness probability Pn defined above and we will 
denote it by the same symbol p„. We may write it as the exponential of an 
entropy difference, 



.AS 



Pn = e- " = e™, (3.2) 

where 5" and 5* are the entropies of the constrained and unconstrained seed 
configuration, respectively. Let ai and 02 denote the area of the annulus and 
the area of the inner disk plus annulus, respectively. As functions of the two 
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lengths R and w they read 



fli = 7r{2Rf -7r{2R-wy, 

as = 7r(2/?)2. (3.3) 
In terms of these two areas we have 

e^ = V^, e'' = (^^^ [L' - a2)^-" a^ (3.4) 
Using fl3.4p in fl3.2p and taking the hmit A^, L — oo as stated we get 

P. = ^e--. (3.5) 

Taking the logarithm of (13. 5p and using (13. 3p we obtain 

logp„ = — log n\ + n log pai — pa2 

= -\ogn\ + nlog[p7iw{4R-w)]- pn{2Rf. (3.6) 

At this point we will take into account the screening effect. Motivated by the 
discussion of section [3?T] we impose that w = 2cRn~'^, where c is a numerical 
constant. This ensures in some average sense that the seeds in the annulus 
are not screening each other. Eliminating w from (13. 6p in favor of R and n 
we find 

logp„ = -logn! + nlog(87rpci?^n~2) -AnpR^ + C(n"^), (3.7) 

valid in the limit n oo. Eq. (13.70 represents the entropy loss associated 
with the constraints on the seed positions. The only freedom that we have 
left is to choose the relation between R and n. Varying (13. 7p with respect to 
R shows that it has a minimum for R = R*, where 

This expression for R*, obtained here heuristically, coincides with the exact 
result of Refs. [4], [Sj for the radius Rc of the n-sided Voronoi cell. The annular 
width w* corresponding to the maximum entropy is obtained by substitution 
of ([31]) in CT, 

w* = c{7ip)~^n~^. (3.9) 

Eqs. (13. 8p and (13.90 show that as n grows, the radius 2R* of the annulus 
of first-neighbor seeds increases but its width w* decreases. The typical 
distance i* between two adjacent first neighbors along this circle is i* = 

11 3 

4:iiR*/n = {4:71 /pn) 2 ~ n~2. From the fact that w* ~ n~2 and hence 
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Figure 7: A Poisson-Voronoi tessellation containing a cell of n = 96 sides. The 
seed density is p = 1. One clearly distinguishes the alignment of the first-neighbor 
seeds along an almost continuous curve. The deviations of this curve from circu- 
larity (see section [5]), predicted to disappear when n — > oo, are still important for 
this value of n. 

<^ <^ 1 it follows that in the large n limit the first order neighbors 
align along an almost continuous curve. This is confirmed by independent 
Monte Carlo simulation according to the method of Ref. [9]. Fig.[7] shows a 
Poisson-Voronoi tessellation containing a 96-sided cell. The alignment of the 
first-neighbor seeds is clearly visible. This figure shows still another feature 
of such a large Voronoi cell: although the cell itself is convex, this need not 
be the case for curve of alignment of the first-neighbor seed positions. 

This is also the suitable place for another side comment. When the 
"seeds" are physical particles having some nonzero diameter a, obviously 
the present considerations cease to be valid when the density of particles be- 
comes too high. From the preceding discussion it is clear that an appropriate 
criterion for their validity is that a -C £*. This point was briefly discussed in 



To complete the calculation of this section we determine the expression for 
the sidedness probability as given by the present approximation. Substitution 
of dXH} in dSZD gives 



Ref. [30]. 



log Pi 



— log n\ + n log 2cn ^ — n + . . . 



—2n log n + n log 2c — ^ 



ilogn-ilog27r + ..., (3.10) 
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in which the dot terms vanish as n ^ oo. 

We now compare Eq. fl3.10p to the exact expansion (11 ■2p . The leading 
order term of (I3.10p is correct and this is the first time that a heuristic 
argument achieves this. We recall that Drouffe and Itzykson [211 132] derived 
the lower bound Pn = —an log n + . . . with a < 2; their conjecture was that 
a = 2, which we now know 01 |2j to be correct [B]. The second term in 
expansion (I3.10p differs from the exact result (11.20 . Indeed, obtaining this 
second term correctly requires a more detailed study of the arrangement of 
the n seeds in the annulus. This will be the subject of section HI 

4 Improved estimate for 

Building on the experience gained above we can now greatly improve upon 
the calculation of p„. In this section will denote the position vector of 
the mth first-neighbor seed and Pm = \Pm\ will stand for its length, with 
m = 1, 2, . . . , n. All exact work on p„ should start from the 2n-fold phase 
space integral 



where x = 1 (o^^ X = 0) if condition for the n seed positions to each 
contribute a perimeter segment is satisfied (or not satisfied); and where A is 
the area from which the remaining seeds must be excluded in order not to 
interfere with this n-sided cell. In this section we will evaluate (14.11) in an 
approximation suggested by the considerations of section [31 The expression 
Xexp(— p^) is a complicated function of the variables of integration that we 
will not need in detail in the present approach. 

4.1 The 2n-fold integral for pn 

We assume that for large n the shape of the n-sided Voronoi cell will be close 
to a circle. Let a new variable P denote twice this circle's radius and hence 
be representative of the linear scale of the cell. We can then split the 2n-fold 
integration in (14.11) into a single one on P and a set of 2n — 1 integrations 
that define the detailed shape and the orientation of the cell. For this set 
we will take the n polar angles of the Pm and the n — 1 independent ratios 
Pm/Pm-i, m = 2, . . . , n of their radial distances. These ratios may in turn 
be expressed in terms of angular variables. Assuming that with negligible 
error we may factorize the integral on P, to be denoted Jrad, and the angular 
ones, we then have 




(4.1) 



Pn I -'rad-'ang ; 

nl 

in which we will consider /rad and Jang separately. 




(4.2) 



12 




All seeds other than the n first neighbors must stay outside the disk of 
radius P, which happens with probability exp(— pvrP^). Hence, extracting a 
scale factor P from all 2n position variables in the integration (14. II) . we have 

/rad = P" P^"-i e-^-^^ = . (4.3) 

Jo 

When this integral is performed, the seed density p disappears from the 
expression for p„, as it had to. 

The integral Jang in (14.21) requires more discussion. It involves 2n — 1 
dimensionless variables of integration for which we can take any set of 2n — 1 
independent angles. Introducing an extra factor (n — 1)! allows us to order the 
n seeds such that their indices correspond to increasing polar angles. A nat- 
ural choice of angles is shown in Fig.[8l it is the set of polar angle differences 
{im^Vm}, where is the angle between two successive mid-point vectors 
Rm-i and Rm , and rj^ the angle between two successive vertex vectors Sm 
and Sm+i- 

Fig. [8] also illustrates how the perimeter segments may be constructed 
successively when the and r/m are given. Two types of construction steps 
alternate. First, given Sm-i and a direction of departure, continue along this 
direction until the segment Sm-iSm spans the angle r]m-i] secondly, when 
arriving at S^, define a new direction of departure by taking an angle of 
to the left. Iterate n times. 
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It is easy to introduce the into (14. ip by passing to a polar coordinate 
representation for the P^; however, the passage from the ratios Pm/Pm-i to 
the T]m leads to the appearance of a nontrivial Jacobian that we will consider 
now. It follows from Fig. [3] and Fig. [HI where we have introduced an auxiliary 
angle Pm, that 

Pm _ \ORm\ _ COs(C,n - /3m) _ ^ - fan /9 (A A) 

Pm-1 \ORm~l\ COSPrn 

Fig. [8] shows that the angles rjm-i and Pm vary with the position of Sm at a 
common rate. To find the Jacobian we therefore calculate the derivative 

i) sin^ 

the last equality being valid asymptotically for small C,m and jSm, which is the 
relevant limit when n gets large. By going around the perimeter one thus 
finds that the transformation from the Pm/Pm-i to the rjm is accompanied 
by a Jacobian C,iC,2 ■ ■ - ^n- This result quantifies our observation at the end 
of section 13. ![ namely that there is an effective repulsion between the mid- 
points. 

The integral on the angular variables now becomes 



-^ang 1) • -^vcrt-^midp 5 (^'6) 



in which 



'vert 



2tt I J}^ \ ^^,^^,n-l 

ri- 1)! ' 



s.m=l 



'•27r /J^ \ i^^\2n-\ 



/midp = / d^i . . . djn (^1 . ■ . ^n) 1^ I - 2 

Vm=l 



-277 1 = ^— ^ — -, (4.8) 
' (2n-l)!' ^ ^ 



where the delta functions enforce obvious sum rules. Due to these constraints 
the angles and rim will typically take values of order rT^ with average 
l-njn. The exact value (here 277) of the upper integration limits in (14. 7p and 
(14. Sp does not affect the outcome of the integrals to the order in n that is 
relevant for our discussion. 



4.2 Improved result 

We now know all the factors that go into the final result. Substituting (14.70 
and (14. 8 p in (14. 6 p and using the result together with (14. 3 p in (14. 2p gives 

1 ( ^^^ (^-1)' (2^)""' (277)^"-^ 

Vn = — : X (n — 1) X ^ — X 7 — X 7 — 

^ n\ ^ ' 277" {n-l)\ (2n-l)! 

(4.9) 



^772)" 



4772 (2n) 
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This improved estimate for p„ is the end result of this section. The same 
expression occurred at an intermediate stage of the calculation in Ref. [5], 
where it went by the name of pn^ . Eq. fl4.9p has the n dependence that we 
know is exact: upon taking the logarithm we now recover the first three terms 
of expansion (11.21) . The exponential (Svr^)" and the factorial denominator in 
fl4.9l) are therefore correct. However, our derivation has been too crude [53] 
that we may trust the proportionality constant (47r^)~^. This constant turns 
out to differ from the correct one appearing in (11.11) by the factor C. In the 
next section we will briefly discuss the origin of this factor. 

5 Elastic deformations and Gabriel neighbors 

In the preceding section we used that the cell perimeter is close to a cir- 
cle. However, since each of the first-neighbor coordinates was integrated 
over a two-dimensional domain, the integral does not respect this circularity 
property exactly. This section deals with the deviations from circularity. It 
will turn out that the quantitative effect of these deviations on the sidedness 
Pn is beyond the heuristic approach of this paper. 

5.1 Random acceleration process 

In this subsection will denote the position vector of the mth mid-point 
and Rm = |R.m| its length, for m = 1, 2, . . . , n. The ratio between two suc- 
cessive mid-point distances can be expressed as Rm/Rm-i = cos7m/ cos/S^ , 
with the auxiliary angles j3m and 7^ defined in Fig. [HI Hence 

Rm-lRm+l COS 

COS 7m+l /r 1 \ 

R"^ cos 7m COS (3m+l ' 

We set 

R^ = R, + SR^ , (5.2) 

where Rc = R* is given by (13. Sp . and expand both sides of Eq. (15. ip using 
that Prn, 1m, and 6Rm/Rc are all small as n gets large. To leading order this 
yields the second-order difference equation 

SRm-l - '^5Rm + 5Rm+l = - ^ Rc {Pm - 1m + 1m+l ~ Pm+l) 

= ~\Rc [Cm (/3m — 1m) — ^m+l{Pm+l — 1m+l)] , 

(5.3) 

where we made use of the identity + 1m = im- In the second line of 
(15.30 the differences (3i — ii (with i = m,m + 1) are of zero average. In that 
line we now replace and ^m+i by their average 27i/n, the idea being that 
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the random contribution to these angles will have a negligible effect when 
Eq. (15.31) is integrated (that is, summed on m) [31]. The right hand side of 
Eq. (15.31) having thus become linear in the angles, it appears that we can 
express the sum of the /3's and 7's again in terms of the ^'s and ?7's. In that 
way we obtain 

i_ 

SRm-l - 2SR^ + 6Rrn+l = (5.4) 

[pn) 2 

with 

The left hand side of Eq. (15.41) represents the "radial acceleration" of the 
perimeter as it turns around the central seed. The right hand side is a random 
term of known properties. An equation like (15.41) also appeared in Refs. [5] 
and [T3], where the conditions of its validity were discussed. It was shown 
there that in the limit n —>■ oo, with the replacement m i— *• = 27rrrm~^, the 
equation becomes a second-order differential equation d^5i?((^)/d0^ = -F(0) 
where F{(f)) is Gaussian noise. For white Gaussian noise this equation is 
referred to as the random acceleration process [IHl EHl [IB] • 

For our present discussion Eq. (15. 4p suffices as it stands. The key point 
is that Fm and F^' are essentially uncorrelated when \m — m'\ > 2. "Es- 
sentially" here and below will mean: neglecting the weak anticorrelation, of 
relative order 0{n~^), induced by the delta function constraint in the inte- 
gral in (14.81) . The F^, have zero mean and are, just hke the and the rj^,, 
of order n^^ as n gets large. 

Let us now sum both sides of (15.41) from m = 1 to some value of m which 
is of order n. The right hand side becomes a sum of 0{n) zero-mean and 
essentially uncorrelated terms that are each of order n~^; hence it will be of 
order n~2. Taking also into account the prefactor {-k / pn)^ we find for the 
"radial speed" 5Rm+i — SRm of the perimeter the expression 

{6Rm+i - 6Rm) - {SRi - 6Ro) = 0{n-'). (5.6) 

Since the average speed must be zero, this means that also 6Rm+i — SRm = 
C(n-i). 

Being the sum of essentially independent increments, the radial speed 
6Rm+i — SRm varies with m as a random walk trajectory. Such a trajectory 
has strong long-range correlations. Therefore, if we sum (15. 6p once more over 
0{n) values of the index m, we obtain the sum of 0{n) strongly correlated 
terms that are each of order n~^, which yields 

6Rm-SRi = 0{l). (5.7) 

We have argued in section 13.21 that the cell perimeter is contained in an 
annulus of width ~ n~2. Eq. (15.71) now shows that this annulus is not strictly 
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Figure 9: Pair of solid lines: the narrow annulus of width w* ~ n~2 through 
which the cell perimeter runs. Its radius has variations of order nP with respect to 
an average radius i?c ~ nz and stays, typically, between the two concentric dashed 
circles. 

circular but that its radius deviates from the typical value Rc = (n/Airp)^ 
by amounts of order unity. In full coherence with this, an evaluation of the 
radial integral (14. 3 p by the steepest descent method yields a peak width also 
of order unity. The resulting large-scale behavior of the cell perimeter has 
been schematized in Fig. [91 A "real-life" example of a not-quite-circular cell 
is provided by Fig. [71 

5.2 Origin of the constant C 

In section [Sm we identified the random acceleration process as the mechanism 
that governs the perimeter's deviations from circularity. These deformations 
may be referred to as "elastic," with the understanding that this elasticity is 
of purely entropic origin. It is in fact a long-range, "macroscopic" elasticity 
that comes over and above the short-range perimeter fluctuations already 
taken into account in section [H by the integrations on the and rjm- 

The exact determination of the effect on p„ of these elastic deformations 
is beyond the heuristic arguments of this paper. For completeness we recall 
that they may be treated [H [3] as a superposition of independent Fourier 
components with wavenumbers q = ±1,±2,.... Their amplitudes ipq are 
Gaussian random variables and their elastic energy [29] E = n Kqipqifj^q 
has the dispersion relation 

Kq = l-q-^ + Aq-\ (5.8) 
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Figure 10: A cell perimeter (heavy solid line) for which the mid-point Rm does 
not lie on the mth perimeter segment SmSm+i but elsewhere on the line prolonging 
this segment. The angle (3m+i is negative. The and are typically of order 
n~2 and the situation depicted in this figure is generic rather than exceptional. 

The elastic deformations multiply p„ by an extra contribution equal to their 
partition function C, which is obtained as C = H^i^g^- 

5.3 Gabriel neighbors of the many-sided cell 

In this subsection we apply the heuristic arguments developed above to 
Gabriel neighborship. Two adjacent Voronoi cells are called Gabriel (or full) 
neighbors if the line segment connecting their seeds does not intersect any 
other cell. In the notation of this work, the condition for the mth first neigh- 
bor of the central seed to be a Gabriel neighbor is that the mid-point Rm 
lie on the perimeter segment SmSm+i- In the example of Fig. [10] the mth 
neighbor is not Gabriel. Reasons for being interested in Gabriel neighbors 
have been mentioned in the Introduction. 

Whereas in two dimensions the average total number of neighboring cells 
is exactly (n) = 6, the average number of full neighbors is equal only to 
^^^fuU _ g^g]^ now: given an ra-sided cell, what is its average number 

u^^^ of Gabriel neighbors? The large-ra behavior of i/^^" can be estimated as 
follows. 

Fig. [To] shows that Rm is in the interval SmSm+i if and only if 7^ and Pm+i 
are both positive. These conditions are not independent since ■jm + Pm+i = 
rjm+i and rjm+i is necessarily positive. An equivalent condition for the mth 
first neighbor to be Gabriel is that 7m > and ?7m > 7m- The advantage of 
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this formulation is that 7m and ?7m are essentially independent (in the same 
sense as before); this is also clear from the algorithm described in section H] 
by which the perimeter is constructed on the basis of given sets of ^'s and 
r^'s. Let us write w„(7m) and Vnijjm) for the probability distributions of the 
7m and the ?7m, respectively. We then have 

/•oo /'OO 

yT = n / d7W„(7) / dr/t;„(r/), (5.9) 
Jo 

in which the prefactor n comes from the total number of cell sides and the 
remaining integral represents the probability that •jm, Pm+i > 0. We now 
need the distributions f„ and Wn- Eq. (14 .yp shows that the rjm define the 
partition of a circle into n random and independent intervals of average 
length 2Tx/n. Hence Vn{f]) is an exponential of average 27r/n, 

v^{ri) = (n/27r)"^exp(-nr//27r). (5.10) 

In order to find Wn we remark that any difference 7^2 ~ 7mi can be written as 
7m2 - 7mi = YZ=mi+ii^rn - Vm-i) , which is a sum of 2m2 - 2mi essentially 
independent random variables. It follows that the differences 7^2 — 7mi , as 
well as the 7m themselves, are Gaussian distributed with a root-mean-square 
width that scales as n~2. Knowing that the 7m must have average vr/n and 
writing the asymptotic behavior of their variance as ~ (con)~^ with cq a 
numerical constant, we have that in the large-n limit 

Wnil) = (con/27r)5 exp [-icon(7 - ir/nf] . (5.11) 

Of course the angles /3m, related to the 7m by inversion of the orientation 
of the perimeter, also have the distribution Wn- We refer again to Fig.fTO] 
for an additional remark. Since Wn is a distribution on the scale n~2 but 
the angles SmOS^+i = Vm are of order n~^, the situation depicted in this 
figure is generic: the perpendicular from the origin onto the line containing a 
perimeter segment will typically hit that line in a point outside that segment. 

Substituting (I5.10p and (15. lip in (15. 9p and evaluating in the limit of large 
n we obtain the scaling with n of the number of Gabriel neighbors, 

z/f"^^ ~ (27rcor2)^ n ^ 00. (5.12) 

This confirms that for many-sided cells Gabriel neighbors are atypical. The 
proportionality constant Cq in (15.120 cannot be obtained by the present 
heuristic arguments; it is, however, in principle accessible by a detailed cal- 
culation in which, again, the Ag will intervene. 

6 Conclusion 

We have developed heuristic arguments applicable to Voronoi cells and whose 
validity is confirmed by existing exact calculations. The arguments, which 
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concern Voronoi cells in the limit of large sidedness, are based on the esti- 
mation of the entropy of a set of points arranged in a plane under specific 
geometrical constraints. The power as well as the limitations of these ar- 
guments have been indicated. As a new result we have determined by the 
reasoning of the present paper the asymptotic n dependence of the average 
number of Gabriel neighbors of the n-sided cell. 

We may henceforth apply these heuristic methods with a certain confi- 
dence in other contexts where no exact calculations are possible. Examples 
of related problems in mathematics and statistical physics have been men- 
tioned in the Introduction. Finally, this work provides the tools that will be 
used in an upcoming study of Voronoi cells in higher spatial dimensions. 
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